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The expected value and the variance of the Baseband Assembly symbol signal-to-noise 
ratio (SNR) estimation algorithm are derived. The SNR algorithm treated here is desig- 
nated as the Split Symbol Moments Estimator (SSME). It consists of averaging the first 
two moments of the integrated half symbols. The SSME is a biased , consistent estimator . 
The SNR degradation factor due to the jitter in the subcarrier demodulation and symbol 
synchronization loops is taken into account. Curves of the expected value of the SNR 
estimator versus the actual SNR are shown. 


I. Introduction 

The Baseband Assembly 1 uses a Split Symbol Moments 
Estimator (SSME) algorithm to estimate the symbol signal- 
to-noise ratio (SNR) of the input signal. Here we describe the 
SSME algorithm and give the expected value and the variance 
of the SNR estimator. Two numerical examples corresponding 
to the Voyager and the Pioneer missions are included to 
illustrate its performance. As in previous Baseband Assembly 
analyses (Refs. 1 and 2), Nyquist sampling rate is assumed. 

II. Statistics of the SNR Estimator 

Figure 1 is a How chart representation of the SSME algo- 
rithm. Referring to this figure, the input to the SNR estimator 
is a string of signal samples modeled as 

y.. = s.. + n.. (1) 

J l] i) i) V A 7 


1 Deep Space Network/ Flight Project Interface Design Handbook , JPL 
internal document 810-5, Rev. D, 1981. 


where 

s.. - unitless random variable whose amplitude is propor- 
tional to the information signal voltage 

n.. ~ unitless random variable whose amplitude is propor- 
tional to the rms noise voltage 

The first two moments of y are 


E{y t ) = s/s 

( 2 ) 

Eiiyy) = S + O" 

( 3 ) 


It is assumed that E {n^ } = 0 

i = 1,2,...,^ Nyquist samples per symbol 
/ = 1,2,...,/? symbols 
The variance of the noise process is assumed to be 
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where N 0 is the one-sided noise spectral density, and B n is the and, consequently, the statistics of Yj will be equal to those of 
one-sided baseband noise -equivalent bandwidth. 


As shown in Fig. 1, in the upper “arm” the samples from 
the first half of a symbol are summed to produce Y oj . In the 
lower “arm” the samples of the second half of a symbol are 
summed to produce Y^. In this analysis, it will be assumed 
that the number of samples in both summers are equal at the 
instants when Y a - and 4 are sampled. For this reason, Y a j 
and Ypj have identical statistics. Making Y aJ - = Y^ = Yj, the 
mean value and the variance of Yj will be, assuming that the 
samples are independent, 


In the SSME algorithm, the random variables Y a and Yp 
are combined to create two new random variables X p and 
X ss in the following way: 


Y Y a 

QL 0 


( 10 ) 


and 


= (»; + V 2 


(id 


Y. = E{Y a .} = E{Y 0j }=^-VsVd. (5) 

- V 2 > - E ^i - = T 


The factor dj , designated as the SNR degradation factor, is 
due to the phase jitter and timing jitter in the subcarrier 
demodulation and symbol synchronization loops, respectively. 
In general, 


0 < dj < 1 


It can be shown that 


d. = 
/ 




(7) 


( 8 ) 


where 

0^. - phase error in the subcarrier demodulation loop 
during the jth symbol 


Tj = timing error in the symbol synchronization loop 
during the /th symbol 


Then, as shown in Fig. 1 , n samples of X p and X ss are averaged 
in the second pair of summers to produce m p and m ' ss . Finally, 
m p and m ss are scaled and combined to produce the random 
variable R*, which is the SNR estimator of the SSME algo- 
rithm, namely, 



The statistics of R* can be determined from the statistics of 
the random variables along the two paths in Fig. 1. These 
statistics are obtained in what follows. 

Using Eqs. (6) and (7) and the fact that Y a and Y@ are 
independent, the first two moments of their product defined 
in Eq. (10) will be 


x p = N 2 Sdl 4 (13) 

X 2 = (N 2 s Sd/4 + No 2 J2) 2 (14) 

The first two moments of X ss defined by Eq. (1 1) are obtained 
using Eq. (A-2) of Appendix A with ^ = N s y/Sd and a 2 = 
A^a 2 , namely, 


X =E{X\) 

ss v s ' 


p T = probability of symbol transition 
T = symbol time 

In this preliminary analysis, it will be assumed that there is no 
doppler stress in the tracking loops and that 0 ; . and r y - are 
functions of the phase and timing jitter only. With this assump- 
tion, (pj and Tj will be constant during one update interval, 
and the subscript j can be dropped, i.e., we will assume that 
during the estimation interval 


= N 2 s Sd + N s o 2 n (15) 

4 = *<4) 

= 3N 2 a 4 + 6N 3 Sd a 2 + N 4 S 2 d 2 (16) 

s n s ns v ' 

Referring to Fig. 1, and using Eq. (A-5), the first and second 
moments at the outputs of the second pair of summers will be 


d , - <*/., 


= d 


(9) 



(17) 
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K) 2 - x l 

(18) 

m 'ss = X ss 

(19) 

K/ = x l 

(20) 


The variances of m ' p and m’ ss are obtained using Eq. (A-6) 
with the moments obtained in Eqs. (13) through (16), namely, 

var(m') = - (X 2 -(X ) 2 ) 
v p 7 n y P v P J } 

N 2 o 2 

-^r ( 21 > 

var K*> = 7 

7N 2 o 2 

= _^L ( 2N s Sd+o 2 n ) (22) 

The covariance of X ss and X p is, by definition, • 

= B{iX„-X2)V p -X p )} (23) 

= E{[ l Y c Y e-r a fy 
= E{[Y 2 a -Yl)HY 2 0 -Y 2 ) 

+ 2(Y Y - TTa] [Y Y.-TTA) 

V a. 0 a 0 /J 1 a: (3 a (3 J J 

+ 2(^^-(r a ) 2 (V 2 ) (24) 


Having obtained the moments of and ra^, we now are 
ready to determine the statistics of the estimator R*. Using 
Eq. (A-9), the expected value of R* defined by Eq. (12) is 


R* = R* 


+ T 


m 


m 


— var(m')+-^ var K»> 


3m' 2 

p 


3 m 


m 


m 


3 2 R* , , , , 

p ss 


m 


m 


(27) 


Inserting Eqs. (A-16), (A-18), (A-19), (21), (22), and (26) in 
Eq. (27) and ignoring higher order terms, we obtain 


R* = R + — (2 R + 1) 
n 


(28) 


where 



is the degraded symbol SNR at the input to the SNR esti- 
mator. From Eq. (28) we observe that R* is a biased but 
consistent estimator (i.e., the bias goes to zero when n goes to 
infinity). 


The variance of R* is obtained using Eq. (A-10), namely,-, 


Using Eqs. (6) and (7) in Eq. (A-2), the third moment of Y a 
and Yq is 

K = ^ = j "I y/Sdol + j A 3 (Sd) 3 ' 2 (25) 

Inserting Eqs. (6), (7), and (25) in Eq. (24) and dividing by 
n , we obtain the covariance of m p and namely, 


var (R*) 



var(Wp) + 



2 


var K*) 


m 


m 


m 


r 

m 


ss 


cov {m. 


m 


N 2 a 2 

s n 

In 


(2NS. + a 2 
x s a n 


(26) 


+ 2 


dm' dm 1 

p SS 


cov (m ' p , m' ss ) 


(30) 
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where 


Inserting Eqs.’ (A-15), (A-17), (A-19), (15), (16), and (26) in 
Eq. (30), we obtain 


var (R*) = (1 + 4R + 2R 1 ) (31) 

By defining the SNR of our estimator as the ratio 

( R *) 2 

SNR (R*) = - 1 > (32) 

var(tf*) 

we see that 

lim SNR (R*) = — (33) 

R-+ o n 

lim SNR (R*) =^ + 2 (34) 


III. Evaluation of d 


B L j - one sided noise-equivalent bandwidth, / = 1 for 
subcarrier loop, / = 2 symbol synch loop 

T l = loop update time, assumed to be identical for both 
loops 

K = number of symbols between updates 
r = symbol time = l/r 

EjN Q = ratio of energy per symbol to noise spectral 
density 


A 

= R = 




(37) 


= M/N s = ratio of the width of the middle portion of 
a symbol to the total symbol length 
(typically 1/2) 

a 2 = L/N ~ ratio of the width of the transition portion 
of a symbol to the total symbol length 
(typically 1/4) 


Assuming that there are no doppler or quantization errors, 
the SNR degradation factor defined in Eq. (8) is a function of 
the phase jitter in the subcarrier demodulation loop and the 
timing jitter in the symbol synchronization loop. Both jitter 
processes, 0 and r, are modeled as Gaussian random variables 
having zero mean and variance and a^, respectively. 


According to Ref, 2, the variance of the phase error in the 
subcarrier demodulation loop at update instants is 


The expected value of d in Eq. (8) will be 


L exp (4 

x v^L ('-^) ! “ p (-ig* v*> 




(35) 



Repeating the steps of Ref. 2, it can be shown that the variance 
of the timing error in the symbol synchronization loop at 
update instants is 



(39) 



(36) 


In Appendix B two numerical examples are given for 
parameter values typical of the Voyager and Pioneer missions. 

In general, the bias in /?*can be reduced by increasing n 
(number of symbols in the estimator). Of course, we can 
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improve our knowledge of R if we compensate for the effects 
of the bias and the degradation factor in Eq. (28), i.e., we may 
assume that the actual input SNR is 



where (R *} is the average value of many R* and d is our esti- 
mate of d. 


IV. Conclusions 

In this article the expected value and the variance of the 
SSME SNR estimator was derived. This estimator was shown 
to be biased and consistent. 

Figures 2 and 3 illustrate the numerical results for the 
Voyager and Pioneer missions. At high signal SNR, the posi- 
tive bias of the estimator dominates over the degradation 
effect due to phase jitter in the tracking loops. At low SNR, 
it is the other way around. Figure 4 is for the ideal case when 
there is no jitter in the tracking loops (d - 1 ). 
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Fig. 3. Mean value and SNR of SNR estimator vs actual SNR: 
Pioneer mission 
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Appendix A 
Gaussian Moments 


1. Relation Between Statistical and 
Probabilistic Moments 

Given a random variable x with Gaussian pdf G(p..a 2 ) and 
defining the rth moment as 

4 = E {*0 (A-l) 


with means 


E {.Kj,} 0 K 

e Ae v e 2 ,...,e K 

it can be shown (Prob 10.17, Ref. 3) that 


(A-8) 


the first four probabilistic moments of x will be 


= M 


t _ 2.2 

M + o 


= 3ff 2 /r + fi 3 


= 3a 4 + 6 o 2 ii 2 + /x 4 


(A-2) 


Defining as the rth statistical moment of a random variable 
, 1 


m 


~n X {x i f 
/'=! 


and the variance of as 


var K) = £ ji 


Chapter 10 of Ref. 3 shows that 

E {m' r } = n' r 
and 

var (m^) = [v 2r - (M' r ) 2 ] 
This is an exact result. 

Given a function g of A' random variables x k , 

#(x) = 


(A-3) 


(A-4) 


(A-5) 


(A-7) 


i * 

E {g (x)} = g(6 ) + Y 53 — r^ (x) 




var (x. ) 


dg(x) 


K K 

H dx.dx. 
1= 1 /= 1 / / 

i*j 


cov (x„ X.) + 


X = 0 


(A-9) 


The variance of £(x) will be (Eq. ( 1 0. 1 2) of Ref. 3) 
K 


var 


= Zj [‘ a ^^ (x) ] 2 I var( ** } 


x - 0 


i= 1 /= 1 1 > 


i * j 


cov (x., x.) + . . . 
v f r 

X = 0 

(A- 10) 


2. Evaluation of the Derivatives 
of the Estimator 

In the SSME algorithm R* is computed from m ' and 
m ss , namely, 


m 


R* = 


2 (-7 m' - m') 
\4 pi 


(A-l 1) 


(A-6) It can be shown that 


and 


E {m' p } = m p = j N 2 Sd 




(A-12) 
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Let 


1 


(A-17) 




E _ K> 


N Sd 

S 


2<, » 


(A- 14) 


The following derivatives of the estimator^* are evaluated: 


a r* 

bm 

p 


a 2 R* 


_ q 1 N 

> > n s 

, m 

p ss 


(1 + 2 R) 


(A-15) 


dm' 2 

P 


(«>,) 


(1+2 R) (A-16) 


m , m 
p ss 


bR* 
a m 

ss 



(A- 18) 


(A-19) 
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Appendix B 

Constants Used to Derive Figs. 2-4 


In order to illustrate the performance of the SSME esti- 
mator, two cases are considered. 

(1) Voyager 

Data rate r = 20,000 symbols/second 

Update time T L = 2.5 seconds 

K = n = 2.5 X 20,000 symbols/loop 
update 

T l B l = 0.4153 (from Table 1, Ref. 1 
for both tracking loops 

Noise bandwidth^ = 3.75 MHz 


(2) Pioneer 

r - 8 symbols/second 
T l = 2.5 seconds 

K = n = 2.5X8 symbols/loop update 
T l B l - 0.4153 for both tracking loops 
5 = 135 kHz 

n 

The performance of R* for the Magellan mission will be 
better than for Voyager. 

Using Eqs. (35), (36), (39), (28), (29), (31), and (32), 
Figs. 2 and 3 are obtained. 
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